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1.  INTRODUCTION 


Digital  filtering  techniques  are  very  appealing  where  complex 
signal  processing  is  desired.  A major  problem  with  digital  filters 
however  is  the  effect  of  quantization  errors  due  to  finite  precision 
arithmetic.  This  thesis  is  concerned  with  roundoff  errors(due  to 
truncation  or  rounding)  occuring  after  each  multiplication  in  linear 
shift-invariant  fixed-point  recursive  digital  filters.  The  accumulation 
of  these  errors  appears  as  noise( error  uncorrelated)  or  limit 
cycles( correlated  errors)  at  the  filter  output  thus  limiting  the  dynamic 
range  of  the  filter  operating  with  a given  wordlength.  This  paper  will 
be  concerned  with  uncorrelated  errors. 

Any  transfer  function  expressible  as  a rational  fraction  in  z can 
be  realized  by  many  filter  configurations.  The  effect  of  roundoff 
errors  however  is  dependent  upon  the  particular  structure  used.  Thus 
structures  are  sought  which  decrease  the  effect  of  roundoff  errors. 
Such  structures  often  lead  to  increased  hardware  costs  and  this  will 
also  be  considered. 

Long[l]  presents  the  design  of  a lowpass  Chebyshev  filter  using  a 
multiple  feedback  structure  which  is  significantly  less  sensitive  to 
roundoff  errors  than  the  popular  cascade  structure.  He  was  unable  to 
get  similar  results  with  a bandpass  design  using  the  multiple  feedback 
structure  however.  In  chapter  3 of  this  thesis,  Long's  lowpass  filter 
is  transformed  to  bandpass  and  to  narrowband  lowpass  in  an  effort  to 
obtain  other  multiple  feedback  filters  with  low  roundoff  noise  output. 
Also  an  equivalent  multiple  feedback  structure  using  less  hardware  is 
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presented.  Chapter  2 discusses  the  modeling  of  roundoff  errors  as 
uncorrelated  random  noise  for  analysis  purposes,  discusses  the  proper 
scaling  of  filter  structures,  and  gives  a direct  relationship  between 
certain  network  sensitivities  and  roundoff  noise.  Chapter  4 analyzes 
the  second  order  coupled  form  structure  which  can  give  better  roundoff 
noise  performance  than  the  popular  direct  form  realization  of  a complex 
pole  pair.  In  chapter  5,  the  numerical  analysis  technique  used  to 
obtain  the  noise  output  of  a filter  structure  is  discussed.  Each  of 
four  transfer  functions  is  realized  by  four  different  structures  and  the 
roundoff  noise  is  determined  for  each. 
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2.  ROUNDOFF  NOISE,  SCALING,  and  SENSITIVITY 
iL  MQdelljK  Roundoff  Noise 

In  order  to  determine  the  roundoff  noise  properties  of  a digital 
filter  structure,  the  roundoff  errors  will  be  modeled  as  random  noise 
sources  injected  at  the  outputs  of  infinite  precision  multipliers  with 
the  following  assumptions: 

1.  The  error  sequence  from  a noise  source  is  a sample  function  of  a 
stationary  random  process. 

2.  Samples  from  the  same  noise  source  are  uncorrelated  with  each  other 
and  are  uncorrelated  with  samples  from  other  noise  sources  and  the 
input  sequence. 

3.  The  probability  distribution  of  the  error  process  is  uniform  over 
the  quantization  range. 

This  method  of  analysis  has  been  shown  to  give  good  results  when  the 
quantization  level  is  not  low  and  when  the  spectral  content  of  the 
filter  input  is  fairly  high[2].  Since  the  filter  is  linear  shift 
invariant,  the  total  output  noise  variance  can  be  determined  by 
superposition. 

Each  noise  source  results  from  the  quantization  of  a multiplier 
output  word  and  its  noise  variance  is  given  by  (2.1)  for  the  uniform 
probability  density  function  where  d is  the  quantization  increment. 
Rounding  will  be  assumed  so  that  the  mean  is  zero. 

o2q  - d2/12 


(2.1) 
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B.  Scaj ipg 

Before  the  noise  variance  of  a filter  can  be  calculated,  overflow 
constraints  must  be  imposed  and  met  at  each  summation  point.  This  is 
accomplished  by  scaling  appropriate  multipliers  in  the  filter.  The 
theoretically  calculated  variance  of  an  improperly  scaled  filter  can  be 
radically  different  from  that  of  the  same  filter  properly  scaled,  for 
example,  consider  the  following  lowpass  filters. 


1-3 

x— O^)- 


3T- 


-i 


L-O—1 


Figure  2.1b.  Filter  b 


The  transfer  functions  of  both  filters  are  given  by  (2.6)  yet  their 
noise  variances  due  to  the  feedback  multiplier  p are  as  given  in  (2.7a) 
and  (2.7b)  corresponding  to  Figures  2.1a  and  2.1b  respectively.  Filter 
b has  been  scaled  in  a manner  to  be  described  below. 


The  noise 
filter  a 


H(z) 


Uzl 

X(z) 


1-3 

1-3 


-1 


o<  3<  l 


a 2 

_a  = xker. 


o 

i 

.2 


1-3 


1-3" 


(2.6) 


(2.7a) 


(2.7b) 


output  of  filter  a is  less  than  that  of  filter  b,  however 
is  improperly  scaled  allowing  overflows  to  occur  at  the 
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summation  node  and  therefore  is  useless. 

Various  methods  for  scaling  filters  exist  and  each  yields  slightly 
different  scaling  multipliers[2,3] • Scaling  in  this  thesis  is 
accomplished  by  satisfying  the  constraint 

max  U (ejU)|  = 1 (2.8) 

-TI  < UU  < TT 

at  each  node  i where  a summation  occurs.  G^i(e^x)  is  the  transfer 
function  from  the  input  to  the  summation  node  i.  This  is  the  L -norm 

CO 

discussed  by  Jackson[3].  Filters  scaled  with  norms  other  than  that  of 
(2.8)  will  have  slightly  different  noise  properties,  however  general 
trends  in  the  noise  performance  of  a structure  should  be  about  the  same. 
Two  s complement  arithmetic  is  assumed  so  that  overflows  may  occur 
internally  at  a summation  node  with  more  than  two  inputs  while  the  net 
sum  still  meets  the  requirement  of  (2.8). 

cu  Jiss  Of  Extra  Scaling  Multipliers  io  Reduce  Noise 

It  was  demonstrated  above  that  scaling  affects  the  output  noise 
variance  of  a filter.  It  is  possible  to  improve  the  noise  performance 
of  some  scaled  filter  structures  by  selectively  introducing  extra 
multipliers  in  the  circuit.  This  will  increase  the  cost  of  the  filter, 
but  if  power-of-two  shifts  are  used  to  accomplish  the  scaling  as 
described  by  Jackson[2]  the  increase  in  cost  will  not  be  as  much  as  it 
would  be  if  full  precision  multipliers  were  used.  Consider  the  internal 
summation  node  i of  a filter  as  illustrated  in  Figure  2.2a.  Such  a node 
would  occur  in  the  common  cascade  filter  structure  for  instance.  The 
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a's  include  proper  scaling  for  node  i.  The  noise  output  due  to  errors 
introduced  at  this  node  is  given  by 


-f  = (m  + n)  -j  . (2.9) 

a * 

o o 

In  Figure  2.2 b,  an  extra  scaling  multiplier  c is  introduced  so  that 
there  are  two  summation  nodes  j and  i.  The  value  c is  chosen  so  that 
when  the  a's  are  divided  by  c,  the  new  overflow  constraint  at  node  j is 
met  exactly.  By  multiplying  by  c at  the  output  of  node  j,  the  overflow 
constraint  at  node  i is  also  met  exactly  as  before.  The  noise  variance 
from  the  m+n+1  multipliers  in  Figure  2.2b  is  given  by 


(m  + 1 + c n) 


(2.10) 


From  (2.9)  and  (2.10)  it  can  be  seen  that  the  introduction  of  the 
multiplier  c will  reduce  the  output  noise  if 


c < . (2.11) 

The  condition  of  2.11  is  often  met.  The  effect  c has  on  reducing 
the  noise  becomes  greater  as  c becomes  smaller.  For  high  Q filters,  c 
could  easily  take  on  values  of  0.1  or  less.  In  such  a case  the  noise 
introduced  by  the  n multipliers  entering  node  j is  insignificant 
compared  to  that  of  the  m+1  multipliers  entering  node  i.  The  noise  from 
the  m+n+1  multipliers  for  this  case  is  approximately 
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(2.12) 


Figure  2.3  show3  a cascade  of  two  second  order  sections  with  the 
multipliers  of  interest  outlined  by  the  dotted  box.  Here,  n=3  and  m=2. 
If  the  second  section  has  high  0,  then  the  improvement  in  noise 
generated  by  the  outlined  multipliers  over  an  equivalent  filter  without 
the  scaling  multiplier  c is  -2.2 dB  as  calculated  below. 


Improvement  “ 10  log  3 -2.2dB 


(2.13) 


The  overall  improvement  in  total  output  noise  of  a filter  ootained  by 
this  technique  will  be  determined  by  the  rest  of  the  filter  and  may  or 
may  not  result  in  a significant  improvement  in  the  total  output  noise. 


JL  A Relationship  Between  fionndoff  Noise 
and  ggcfcaln  Nstygrk  Sensitivities 

This  section  shows  how  the  roundoff  noise  may  be  expressed  exactly 
in  terms  of  certain  network  sensitivities.  Fettweis[4]  has  developed  a 
similar  relationship  however  it  does  not  express  the  total  roundoff 
noise  exactly. 

Fettweis  has  shown[5]  that  the  sensitivity  of  the  input-output 
transfer  function  Hn1  witn  respect  to  a multiplier  a directed  from  node 
k to  node  i is  given  by 
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dH  1 

-a“-HklHni  <2*“> 

Hni  is  the  transfer  function  from  node  i to  the  output  and  is  the 
function  required  for  the  computation  of  the  gain  of  the  noise  power 
from  node  i.  Solving  (2.14)  for  Hni  and  substituting  this  in  (2.3) 
yields 


o 


II 

ij  l-J- 

2"  J_lHkl 


2 

djL 


(2.15) 


This  expression  gives  the  roundoff  noise  variance  in  terms  of  the 
sensitivity  of  Hn-|  with  respect  to  ot  divided  by  the  gain  Hk1  from  the 
input  to  the  multiplier  source  node  k.  This  gain  Hk1  is  dependent  on 
the  filter  structure.  However  if  k=n,  then  (2.14)  can  be  written  as 

2 

1 ‘ IjT"  'sSIl  da)  (2.16) 

°o  nl 

Hn1  is  the  input-output  transfer  function  and  is  not  related  to  the 
filter  structure.  Thus  the  noise  variance  due  to  a noise  source  at  node 


i is  expressed  exactly  in  terms  of  the  sensitivity  of  a feedback 
multiplier  a from  the  output  to  node  i.  This  multiplier  is  often 
present  in  a structure  although  it  need  not  be  for  (2.16)  to  hold. 
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3.  THE  MULTIPLE  FEEDBACK  STRUCTURE 


At  Motivation 

Tne  multiple  feedback(MFB)  or  leapfrog  digital  filter  structure  is 
shown  in  Figure  3-1.  The  transfer  functions  T^Cz)  are  of  the  form 


Tt(z) 


a0i+aUZ'1 +a2iZ'2 
1-0liZ‘1-02iz"2 


(3.1) 


and  are  realized  by  first  or  second  order  sections.  Long[1]  used  the 
technique  of  numerical  coefficient  matching  in  the  z-domain  by  solving 
simultaneous  nonlinear  equations  to  design  several  filters  of  this  type. 
He  determined  the  noise  variances  of  these  filters  and  compared  them  to 
those  of  filters  with  the  same  transfer  functions  but  realized  with 
different  structures.  In  particular,  comparisons  were  made  to  the 
cascade  of  direct  form  second  order  sections.  He  designed  three  fifth 
order  lowpass  Chebyshev  filters  using  first  order  sections  for  the 

T,(z)'s  as  in  Figure  3.2.  The  a and  o multipliers  are  missing  in 

order  to  make  the  filter  physically  realizable.  The  transfer  functions 
of  these  filters  had  0.5  dB  ripple  and  cutoff  frequencies  of  j , , 

TT 

and  g and  are  designated  LP2,  LP4,  and  LP3  respectively.  For  the 
wideband  LP2  transfer  function,  the  MFB  design  was  about  4 dB  worse  than 
the  cascade  LP2  design.  The  LP4  MFB  and  cascade  realizations  showed  no 
appreciable  difference  but  the  LP8  MFB  realization  proved  to  be  10  dP 
better  than  the  LP8  cascade  structure  indicating  that  for  narrowband 
lowpass  filters  the  MFB  structure  can  significantly  improve  the  signal 
to  noise  ratio  of  a filter.  Long  also  designed  6th  oraer  bandpass 


Figure  3.5  A Modification  of  the  Miltiple  Feedback  Structure 


;««M«wy.ya6ge 
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( 3 • 3 b ) respectively  where  N is  the  order  of  the  filter. 


H (Z) 
a 


-k  ? "i 
z n — 

i-l  Di 


where 


k = N/2 


N even 


= (N-l)/2  N odd 


Hb(z)  = 


rr  _ 
i=l  D, 


The  denominators  of  (3. 3a)  and  (3.3b)  are  given  by 

A = 1-(sum  of  all  individual  loop  gains)+(sum  of  gain  products 
of  all  possible  combinations  of  two  non-touching  loops) 

-(sum  of  the  gain  products  of  all  possible  combinations 

of  3 non-touching  loops)*.. . (3.4) 

The  gain  products  in  (3.4)  for  Figure  3*5a  consist  of  terms  of  the  form 


-1  1 

i’i+lZ  DiDi+l 


for  i odd 


-1  i+1 

0t  D.+1 


for  i even. 


The  gain  products  in  (3.4)  for  Figure  3 -5b  consist  of  terms  of  the  form 


. -1  i 

‘•1+1  DiDi+l 


for  all  i. 


Thus  the  denominators  of  Ha(z)  and  Hfe(z)  will  be  the  same  when  Ni=NW1 
for  all  even  i.  This  is  the  case  for  the  Chebyshev  filter  of  figure  3.2 
since  all  zeroes  lie  at  z=-1.  The  LP8  filter  was  designed  using  the  new 
structure  of  Figure  34.  The  noise  variances  were  computed(as  in 
Chapter  5)  for  both  old  and  new  structures  and  differed  by  only  0.1  dB. 


19 


£.  Freg^engi  IrangformaUgps 

In  order  to  transform  from  a lowpass  filter  of  cutoff  frequency  0£ 
to  a bandpass  filter  centered  at  w with  bandwidth  B the  following 
substitution  is  made  for  z-1  in  the  lowpass  prototyped] 


where 


-1 

z 


-2  -l.k 

z - az  + b 

. -2  -1  . 
bz  - az  + 1 


2<*k 

k+1 


b 


k-1 

k+1 


(3.6) 


and 


cos  w 

c 

cos (B/2) 


k * cot  (B/2)tan(0  /2). 

c 


This  substitution  transforms  the  general  first  order  section  shown  in 

Figure  3 da  to  the  second  order  section  shown  in  Figure  3.6b.  To 

transform  the  MFB  filters  in  Figures  3.2  and  3d  to  bandpass  filters, 

the  substitution  depicted  in  Figures  3 -6a  and  3.6b  is  made  for  each 

first  order  section  occuring  in  the  prototype  MFB  filters.  For  the 

transformed  filter  to  have  no  delay-free  loops,  it  is  required  that 

B = 0 (k=1  ,a=  <*,b=0) . Performing  this  transformation  on  the  filter  of 

c 

Figure  3.2  is  straightforward  and  results  in  the  filter  shown  in  Figure 
3.7  where  use  has  been  made  of  the  fact  that  all  zeroes  are  at  z=-1  in 
Long's  LP8  filter.  It  should  be  noted  that  four  more  full  precision 
multipliers,  two  more  delays,  and  three  more  summations  are  necessary 
here  than  would  be  necessary  for  a 10th  order  MFB  bandpass  filter 
designed  like  the  6th  order  filter  in  Figure  3.3. 


Performing  this  transformation  on  the  modified  MFB  filter  in  Figure 
3.<<  is  slightly  more  complicated.  Each  feedback  path  is  considered  to 
be  a multiplier^  in  Figure  3.6a  and  results  in  two  feedback  paths. 
The  transformed  filter  is  shown  in  Figure  3-d  and  has  the  same  magnitude 
characteristic  as  the  filter  of  Figure  3.7.  This  filter  requires  the 
same  number  of  full  precision  multipliers  as  the  other  transformed 
bandpass  filter(Figure  3.7)  however  it  requires  fewer  delays  and 
summations.  In  an  analysis  of  the  noise  variance  of  a bandpass  filter 
design,  the  noise  variance  of  the  filter  in  Figure  3-8  was  only  0.7  dB 
worse  than  that  of  the  filter  in  Figure  3*7.  Since  the  noise  properties 
of  both  the  old  and  modified  MFB  structures  appear  to  be  similar,  only 
the  modified  structure  and  its  transforms  will  be  considered  in  the 
noise  analysis  of  Chapter  5. 

In  order  to  transform  from  a lowpass  filter  of  cutoff  frequency  & 

c 

to  another  lowpass  filter  of  cutoff  frequency  ou  , the  following 

c 

substitution  is  made  for  z_1  in  the  lowpass  prototyped] 


where 


(3.8) 


This  substitution  transforms  the  general  first  order  section  of  Figure 
3.6a  to  that  of  3.6c.  Using  this  transformation  on  the  modified  MFB 
structure(Figure  3.M)  yields  the  filter  shown  in  Figure  3*9.  It  should 
be  noted  that  this  structure  is  physically  unrealizable  as  it  contains 
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delay  free  loops.  It  is  included  for  comparison  purposes  only  in 
Chapter  5. 
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3.  THE  COUPLED  FORM 


Ai  Structure 

The  coupled  form  realization  of  a complex  pole  pair  as  proposed  by 
Gold  and  Radar[7]  is  illustrated  in  Figure  4 . 1 . 


The  noise  properties  of  this  structure  will  be  compared  to  those  of  the 
direct  form  implementation  shown  in  Figure  4.2. 


°1 


Figure  4.2.  The  Direct  Form 


The  transfer  function  for  Figure  4.1  is 


a a z 

0 l 

H(z)  - ^ 77  (A.l) 

1 - (P1-^2)z  + 01*2  ‘ 


and  that  of  Figure  4.2  is 


1 - 2rcos9z  + r z 


Thus  these  two  filters  realize  the  same  pole  pair  if 


2rcos9 


a and  a are  set  by  the  scaling  required  at  nodes  n,  and  np  in  Figure 


4.1.  There  is  one  free  parameter  available  in  designing  this  filter 


It  is  desired  to  select  that  parameter  to  reduce  the  output  noise 


is  defined  by 


The  sensitivity  function  s' 


The  sensitivity  functions  of  (4.1)  are  found  to  be 
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(4.5d) 


where  z has  been  replaced  by  eJ  to  point  out  the  dependence  on 
frequency.  (4.5c)  can  be  minimized  with  respect  to  3^  and  3^  to  obtain 


- &2  ■ rcosB. 


(A. 6) 


Also,  since  (4.5a)  and  (4.5b)  together  are  symmetric  in  8^  and  89  , this 
seems  a good  choice  for  them.  This  selection  is  given  by  Oppenheim  and 
Schafer [6] . 


C.  Noise  Analysis 


To  perform  a closed  form  noise  analysis,  the  transfer  functions 
nodes  n1  and  n2  to  the  output  are  required.  These  are  given  by 


n.  . _ Q -1  , 2 -2 

1 1 - 2rco80z  +r  z 


(1  -B^'Sz'1 


n2  1 - 2rcoa9z  ^ +r2z  2 


The  transfer  function  from  node  n,  to  the  output  of  the  direct  form  is 


nl  direct  i . 2rcosez'1 +r2z'2 


It  can  be  seen  that  the  magnitudes  of  both  (4.7a)  and  (4.7b)  can  be  less 
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than  the  magnitude  of  ( U . 7c ) . Substituting  (4.7)  into  (2.4)  and  using 
Jury's  tables  of  that  integral[8],  the  output  noise  variances  for  a 
single  noise  source  at  each  node  are  found  to  be 


a 2 2 

"l  a\a+r  ) 

(1  - r2)(l +r4  - 2r2cos26) 

a2  2 2 

(l+P^)(l+r  ) - 43^rcos9 

~2  2 4 2 ‘ 

(1  -r  ) (1  + r -2r  cos20) 


(4.8a) 


(4.8b) 


The  output  noise  variance  for  a noise  source  at  node  n.|  of  the  direct 
form  version  is  given  by 


1 + r 


(1  - r2)(l  +r4  - 2r2cos20) 


(4*9) 


Therefore  (4.8a)  and  (4.8b)  can  be  rewritten  a3 
o2 

ni  2 
■ =0  N 

■J  12 


(4.10a) 


I1-  d+e2  . 4rco|0  e (4.10b) 

°o  1+r 

a2  a2 

4rcos0  then  "2  <N  If  these 
a*  £ 


Thus  if  o1<1,  then  -yl<N2  and  If  3^ 

^ 0 

conditions  are  met,  then  the  noise  gain  from  each  node  of  the  coupled 


form  is  less  than  that  of  the  noise  gain  N2  of  the  direct  form.  The 


total  noise  output  for  each  form  will  depend  on  the  number  of  noise 
sources  at  each  node  in  the  filters.  For  the  coupled  form,  there  are 
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three  noise  sources  at  node  n1  due  to  0^,  and  t and  there  are  two 
noise  sources  at  node  n2  due  to  and  9^ . For  the  direct  form  there 
are  three  noise  sources  at  node  n^  The  total  noise  outputs  of  the 
direct  and  coupled  forms  are  given  by  (4.11a)  and  (4.11b)  respectively. 


coupled 


3c^N2  + 2(1  +32 


4rcos6 
1 +r2 


VN2 


3N 


2 


(4.11a) 


(4.11b) 


Therefore  the  coupled  form  noise  is  less  than  the  direct  form  noise  if 


+ f (i+e2i  - 


4rcos9 
1 +r2 


V 


< 1. 


(4.12) 


Closed  Form  Noise  Calculations 

The  objective  here  is  to  select  the  parameters  3^,  3^ , and  t in 
Figure  4.1  so  as  to  minimize  the  noise  variance  given  by  (4.11a).  Since 
N2  is  a function  only  of  the  pole  locations  (4.9),  it  is  sufficient  to 
minimize  the  following  expression 


NUM  = 3o2  +2(1  +32 


4rcos8 
1 +r2 


*!>• 


(4.13) 


To  do  this,  an  expression  for 
constraints  which  must  be  met  at  the 


must  be  obtained.  The  overflow 
summation  nodes  n^  and  n2  are  given 


» 


wa 


by 
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max  I H (e“ 

-n < uu < n wlx 


V'-V'1*’ 


1 „ -Iju  2 -21w' 

-n<uxn  l-2rcos6e  J +r  e J 


1 (4 . 14a ) 


i a a e*^X 

max  I",  x(eJi)l  = "ax  I P"J— ja, T^l=1  <4-14b> 

-tt  < au  < n 2 -n<aj<rn  1 - 2rcos0e  J +re  ■* 


A closed  form  solution  for  a and  a from  (4.14)  is  not  readily 

0 1 

obtainable.  It  can  be  seen  that  and  will  be  functions  of  0^  and 

not  of  t.  If  a closed  form  solution  for  or  were  available,  it  could  be 

used  in  (4.13)  with  0^  replaced  by  2rcos0  - 0^  and  a closed  form 

expression  for  the  minimum  noise  could  be  found  by  taking  the  derivative 

of  (4.13)  with  respect  to  0^  and  substituting  the  resulting  0^  into 

(4.11a).  Since  these  expressions  were  not  available,  the  following 

algorithm  was  used  to  find  the  minimum  variance  on  a digital  computer. 

Step  1:  Select  a value  for  0^  (02  “ 2rcos0  - 0, ) 

Step  2:  Solve  (4.14a)  iteratively  for  Oq 
Step  3:  Solve  (4.14b)  iteratively  for  O' 

Step  4:  Solve  (4.13)  for  NUM 

Step  5:  Repeat  steps  1-4  until  NUM  is  minimized 
An  alternative  to  selecting  0^  to  minimize  the  variance  of  the 
coupled  form  is  to  set0^  to  unity.  This  eliminates  the  0^  multiplier 
and  reduces  the  number  of  noise  sources  at  node  n1  from  3 to  2.  This 

will  be  referred  to  as  coupled  form  A and  the  result  of  the  previous 

design  algorithm( variance  minimized)  will  be  referred  to  as  coupled  form 


B.  The  variance  for  coupled  form  A is  given  by 
2 a2  a2 

JL  m 2 -li  + 2 — 

_2  2 2 

a o o 

0 0 0 


(4.15) 


f 
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By  selecting  0^  = 1,  the  parameters  @2  and  o^t  are  determined  by  (4.3). 
c»q  and  & are  found  as  in  steps  2 and  3 of  the  previous  algorithm.  Then 
the  variance  can  be  calculated  from  (4.15). 

These  calculations  were  made  for  various  pole  locations  and  the 
results  are  Dlotted  in  Figures  4.3,  4.4,  and  4.5.  Each  figure  is  a plot 
of  the  noise  variance  in  dB  verses  the  angle  0 of  the  pole  location  for 
a fixed  pole  radius  r.  Also  included  is  the  variance  of  the  direct  form 
realization  of  the  filter.  The  plots  go  to  0 =90°  and  they  are 

symmetric  from  90°  to  13O°(01=-1  in  form  A for  0>9O°).  From  the  graphs 
it  can  be  seen  that  for  large  pole  radius  r and  small  pole  angle  9,  both 
coupled  forms  A and  B offer  significant  improvement  over  the  direct 
form.  At  0=0,  coupled  form  B is  36*4  dB  better  than  the  direct  form  for 
r=.99  and  it  is  27.2  and  16.5  dB  better  for  r=.95  and  r=.90 
respectively.  For  small  angles,  coupled  form  B is  seen  to  be  slightly 
better  than  coupled  form  A.  However  for  most  angles  less  than  50°,  form 

A which  has  fewer  noise  sources  is  seen  to  be  slightly  better.  These 

results  are  tabulated  in  Table  4.1.  The  structures  are  listed  in 

descending  order  according  to  their  performance.  0^  is  shown  in 

Figures  4. 3-4. 5. 


Angle  Range 


Rank 

0 

1 

CD 

O 

eo  - 50° 

50°  - 90° 

1 

Form  B 

Form  A 

Direct 

2 

Form  A 

Form  B 

Form  B 

3 

Direct 

Direct 

Form  A 

Table  4.1.  Noise  Performance  of  Second  Order  Structures 


Variance  versus  9 for  Second  Order  Sections  (radius  = .99) 


Figure  4.5  Noise  Variance  versus  9 for  Second  Order  Sections  (radius 
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A plot  of  verses  0^  for  r=.95  is  shown  in  Figure  4.6.  It  can  be 

seen  that  at  0 =90°,  0^=0  and  0^=0  and  coupled  form  8 reduces  to  the 

direct  form.  The  discrepancy  in  the  noise  variances  of  Figure  4.4  at 

this  point  is  due  to  the  fact  that  the  noise  calculation  for  the  minimum 

variance  in  form  B was  made  with  5 noise  sources  while  that  of  the 

direct  form  was  made  with  3-  This  difference  in  noise  sources  also 

explains  why  coupled  form  B did  not  perform  as  well  as  the  direct  form 

for  larger  angles.  It  can  also  be  seen  in  Figure  4.6  that  ft  ft  for 

1 *2 

coupled  form  B which  is  the  choice  indicated  by  the  sensitivity  analysis 
of  Section  B.  Experimentally,  it  was  found  that  if  the  noise  gains  at 
nodes  n1  and  n2  are  weighted  equally  as  in  (4.15),  then  the  minimum 
variance  occured  at  exactly  0^=02=rcos0  . 

The  result  that  coupled  form  A performs  better  than  coupled  form  B 
for  angles  between  0^  and  50°  and  is  not  much  worse  for  angles  less  than 
0q  makes  this  form  appealing  since  it  requires  only  two  full  precision 
multipliers  as  compared  to  three  required  for  coupled  form  B.  Thus  the 
low  noise  performance  of  the  coupled  form  can  be  obtained  with  the  same 
number  of  full  precision  multipliers  that  are  required  in  the  direct 
form. 
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5.  NUMERICAL  NOISE  ANALYSIS 


A.-  Analysis  Technique 

The  analysis  method  used  for  computing  roundoff  noise  is  based  on 
the  signal  flow  graph  representation  of  a digital  filter.  The  actual 
noise  analysis  program  used(NOISE)  was  written  in  FORTRAN  by  Long[1] 
based  on  the  method  of  Jackson[2].  This  method  is  briefly  described 
below. 

Using  the  signal  flow  graph  representation,  any  digital  filter 
structure  with  N nodes  can  be  put  in  the  matrix  form 

Y(z)  = z-1HdY(z)+HcY(z)+U(z)  (5.1) 

where  Y(z)  = Nxl  vector  of  node  output  values 

U(z)  = Nxl  vector  of  node  input  values 

Hc  = NxN  matrix  of  coefficients  for 
branches  with  no  delay 

Hd  = NxN  matrix  of  coefficients  for 
branches  with  one  delay 

The  (m.n)^  element  cf  Hc(or  Hd)  is  the  coefficient  of  the  multiplier  in 
the  branch  with  no  delay(or  with  one  delay)  directed  form  node  n to  node 
m.  (5.1)  can  be  rewritten  as 

r 

(I-z‘1Hd-Hc)Y(z)  = U(z)  (5.2) 

The  complex  gain  from  node  i to  node  k at  any  frequency  uuQ=  uj/T  where  T 
is  the  sampling  rate  can  be  found  by  solving  the  linear  simultaneous 


# 
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equations  (5.2)  with  the  following  substitutions 

(5.3a) 
(5.3b) 


U(z)  = [o. 


.010. 

T 


.or 


and 


ith  element 


z = e 


The  gain  from  node  i to  node  k is  then  given  by  the  kth  element  of  the 

n 

solution  vector  Y(e  ).  A modified  version  of  the  IBM  scientific 
subroutine  SIMQ  was  used  to  solve  the  set  of  complex  linear  equations 
(5.2).  It  was  found  that  for  filters  with  very  narrow  bandwidths, 
double  precision  arithmetic  was  required  in  computing  (5.2). 

To  compute  the  noise  power  output  of  a filter  based  on  the  model  of 
Section  2. A,  (2.3)  and  (2.5)  are  used.  Substituting  (2.3)  into  (2.5) 
and  summing  over  the  network  nodes  rather  than  the  noise  sources  yields 

°2  n 

-7  = 5^1  S(a)dua  (5‘4) 

ao  "n 

where 

N 2 

N0*>)  “ iS1kJHni(ejU)>'  (5'5) 


and  N = number  of  nodes 

kj^  = number  of  error  sources  at  node  i 

= transfer  response  from  node  i to  the  output 

N(iu ) is  the  power  spectral  density  of  the  noise.  All  are 

obtained  at  once  by  substituting 
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u(z)  = [o oir 


(5.6) 


in  (5.2).  By  the  interreciprocity  property,  Hni(e^)  is  given  by  the 
ith  element  of  Y(e^).  The  integral  of  (5.4)  is  approximated  by  the 
trapezoidal  rule  as 


N -1 


(5.7) 


where 


-f-^-[&(0)+N(TT)]  N(^) 


\ = 1 • fj-  ’ 1 = 0,1, ...,N  • 

n ‘ 


A value  of  Np=50  was  used  for  all  examples  listed  in  Section  5.B.  This 
value  for  Np  gave  results  accurate  to  at  least  six  decimal  places  when 
NOISE  was  run  on  second  order  structures  for  which  the  closed  form 
variance  could  easily  be  computed  by  hand. 

A measure  of  the  noise  power  over  a band  of  frequencies  can  be  made 
by  changing  the  limits  of  the  integral  of  (5.4)  to  obtain 

.2 


where 


uu 

aPB  1 ' 

2 ~ “ R *T  »<*>d“ 

°0  “l 


B “ ui  -u^ 


(5.8) 


a>2  ■ upper  frequency  (normalized  w.r.t.  T) 

■ lower  frequency  (normalized  w.r.t.  T) . 

The  integral  of  (5.8)  has  been  normalized  with  respect  to  the  bandwidth 
B and  is  approximated  by  (5.9)  for  computer  evaluation. 


-p'sr  +»<»2>’  ♦ r & »<v 


N -1 

1 t 


(5.9) 


in 


where 


0,1, ... ,N  . 

P 


JL  Exain£i&§  £PI  CfimjESllSon  of  filt£r  Structures 

The  filter  structures  which  will  be  compared  in  this  chapter  are 
the  standard  cascade  structure,  the  modified  multiple  feedback 

structure,  and  the  coupled  form  second  order  sections  in  cascade  form. 
The  coupled  forms  of  Chapter  4 were  modified  by  placing  a zero  after 
each  first  order  section.  This  and  the  fact  that  the  total  noise 
variance  of  a section  in  a cascade  structure  depends  on  the  other 
sections  means  that  the  use  of  coupled  form  no  longer 

guarantees  minimized  noise  power.  It  is  expected  that  the  performance 
will  be  comparable  however. 

Two  lowpass  Chebyshev  filters  and  two  bandpass  Chebyshev  filters 

will  be  used  for  the  comparisons.  All  have  0.5  dB  passband  ripple. 

Also,  all  have  fairly  narrow  bandwidths  since  Long's  results  and  those 
of  Chapter  4 indicate  that  the  MFB  and  coupled  form  structures  will  be 
of  greatest  interest  for  the  narrow  bandwidth  case.  The  magnitudes  of 
the  transfer  functions  are  pictured  in  Figures  5. 1-5.4.  The  two  lowpass 

TT  TT 

filters  have  0.5  dB  cuttoff  frequencies  of  g and  gg  and  will  be 
designated  LP8  and  LP80  respectively.  The  two  bandpass  filters  have  0.5 
dB  bandwidths  of  g . Their  center  frequencies  are  .437"  and  g and 
these  filters  will  be  designated  BP1  and  BP2  respectively. 


The  filters  and  their  coefficients  are  given  in  Figures  5-5-5.11 
and  Tables  5- 1-5.6.  The  coefficients  for  the  MFB  structures  come  from 
Long's  MFB  filter  and  transformations  of  it.  The  ordering  of  the 
cascaded  second  order  sections  in  Figures  5.5,  5.8,  5.9,  and  5.11  is 
based  on  Jackson[2].  All  filters  satisfy  the  scaling  requirement  of 
(2.8).  This  was  accomplished  on  the  computer  using  the  techniques  of 
Section  5. A to  determine  the  maximum  gain  at  each  summation  node.  Also, 
the  technique  of  Section  2.C  was  used  to  reduce  the  noise  of  each  filter 
where  possible.  The  BP2  filter  was  not  constructed  with  coupled  form  B 
since  its  noise  properties  are  similar  to  those  of  coupled  form  A. 

Each  multiplier  in  a filter  with  coefficient  other  than  one  is 
counted  as  a noise  source  including  those  scaling  multipliers  which  may 
not  be  full  precision  multipliers.  The  multipliers  whose  coefficients 
are  one  may  be  replaced  by  a direct  path  between  nodes  in  the  filter. 


Figure  5.1  LP8  Magnitude  Characteristic 
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Figure  5.2  LP80  Magnitude  Characteristic 
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l£ a Filter 


Section  i 

1 

2 

®li 

.86555015 

1 .72815082 

02i 

- 

-.79460336 

a . 
Oi 

.5 

.25 

a 

li 

.5 

.5 

a 

2i 

- 

.25 

Si 

.13444980 

.066452714 

LPdO  Filter 

Section  i 

1 

2 

Bu 

.98587045 

1 .97651463 

B21 

- 

-.97724146 

°0i 

.5 

.25 

a 

li 

.5 

.25 

“21 

- 

.25 

Si 

.014129550 

.000726930 

3 

1 .72669227 
-.91793529 
.25881416 
.51762833 
.25881416 
.14586934 


3 

1 .98965721 
-.99124730 
.25 
.5 
• 25 

.001590090 


Table  5.1. 


Cascade  Realization  of  Fifth-Order  Chebyshev 
Lowpass  Filters  (Figure  5.5) 


(not  physically  realizable) 
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LI 8 Filter 

(Figure  5.6) 
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.73984200 
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.5 

.51741791 
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.505801 16 

.50142226 

.5 

.51741791 

Si 

.27256869 

.21088026 

.20847000 

.25423218 

C 1 , i+1 

-.27027892 

-.10586551 

.18621970 

-.26794146 

L£80  filter  Ifigure  iJi 
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li 
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.098877889 

-.14450219 

"ati,i+l 

.13820345 

.047837684 

.081108833 

.11853412 

Table  5.2. 

Multiple  Feedback 

Realization 

of  Fifth-Order 

Lowpass  Filters 


5 

.73984366 

.52696216 

.52696216 

.40029000 


5 

.97090530 

.5 

.5 

.047181503 


Chebyshev 
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m Filler 


Section  i 

1 

2 

3 

4 
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.76137695 
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.84242903 
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- 
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- 
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- 

Table  5.3. 

Coupled  Form  A Realization  of  Fifth-Order 
Chebyshev  Lowpass  Filters  (Figure  5.8) 

m mat 
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Table  5.4. 

Coupled  Porn  B Realization  of 
Lowpass  Filters  (Figure  5.8) 

Fifth-Order 

Chebyshev 
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BPI  £Ui£r 


i.  Section  i 
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.37311002 
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.60248480 
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-.90956961 

-.84101291 

-2.50940195 

Si 

.065995086 

.053728000 

.21435000 

.28722000 

.10346000 

Table  5.5.  Cascade  Realization  of  Tenth-Order  Chebyshev 
Bandpass  Filters  (Figure  5.9) 
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m Filter 
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2 

3 

4 

5 

8U 

.34796840 

.4 

.4 

.4 

.34796873 

S2l 

-.73984200 

-1 

-1 

-1 

-.73984366 

“Oi 

.50680116 

.50141044 

.5 

.52830427 

.52679691 

a2i 

-.50580116 

-.50141044 

-.5 

-.52830427 

-.52679691 

Si 

.27256869 

.21088026 

.20422700 

.25423219 

.40019600 

■*1,1+1 

.27027892 

.10806750 

.18621970 

.26248184 

- 

atl,i+l 

-.054055783 

-.021613500 

-.037243940 

-.052496368 

* 

212  Filtfir 

Section  1 

1 

1 

2 

3 

4 

5 

eu 

1 .63889532 

1 .88395880 

1.88395880 

1.88395880 

1 .63889690 

62i 

-.73984200 

-1 

-1 

-1 

-.73984366 

a0i 

2.49279940 

2.08444795 

1.541827 28 

2.52734089 

2.50938850 

“21 

-2. 49279940 

-2.08444795 

-1.54182728 

-2.52734089 

-2.50938850 

Si 

.054918343 

.051084467 

.067606500 

.052049375 

.084058000 

”*■1,1+1 

.22638713 

.078527574 

.29496800 

.26122411 

- 

ati,l+l 

-.21325201 

-.073971351 

-.27785377 

-.24606773 

— 

Table  5.6. 

Multiple  Feedback  Realization  of  Tenth-Order 
Chebyshev  Bandpass  Filters  (Figure  5.10) 
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BP1  Filter 


Section  i 

1 

2 

3 

4 

ei 

1 

.11933007 

1 

.015981910 

P 

o 

H* 

.64514000 

. 7913^000 

.67286000 

.74773000 

*11 

.64514000 

-.7913^000 

.67286000 

-.74773000 

Si 

.13390000 

.98267000 

.20547000 

.99777000 

*"1,1+1 

-1.17707805 

- 

-1.30010299 

- 

BP1  Filter (Continued) 

Section  i 

6 

7 

8 

9 

0i 

.24824219 

1 

-.02301361 

1 

*01 

.79164000 

.70407000 

.71017000 

.62546000 

a 

11 

-.79164000 

.70407000 

-.71017000 

.62546000 

Si 

.98014000 

.20438000 

1.00041000 

.71182000 

ti,i+l 

- 

-1.39056036 

- 

-.91981667 

Table  5.7.  Coupled  Form  A Realization  of  Tenth-Order 
Chebyahev  Bandpass  Filter  (Figure  5.11) 
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5 

1 

.64509000 

.64509000 

.32181000 

-1.02131797 


10 

•35779676 
.79302870 
-.79302870 
1 .0463000 
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BP2  Filtfil 


Section  i 

1 

2 

3 

4 

5 

8 . 
i 

1 

.74913762 

1 

.77303993 

1 

a 

Oi 

.51518560 

2.0900057 

.52228205 

1.7309708 

.52198295 

a 

li 

.51518560 

-2.0900057 

.52228205 

-1.7309708 

.52198295 

Si 

.13646288 

.48308654 

.096567230 

.55629730 

.38955737 

ti,i+l 

-.43760739 

- 

-.56913959 

- 

-.18840886 

m Filter (Continued) 


Section  i 

6 

7 

8 

9 

10 

.81142998 

1 

.89883662 

1 

.95943788 

a 

Oi 

1.74215060 

.52175786 

1.7562198 

.52144475 

1.76660582 

"li 

-1.74215060 

.52175786 

-1.7562198 

.52144475 

-1.76660582 

Si 

.55030243 

.479975231 

.55000000 

.50737774 

.55474390 

ti,i+l 

- 

-.094859640 

- 

-.064832620 

- 

Table  5.8. 

Coupled  Form  A 

Realization  of 

Tenth-Order 

Chebyshev  Bandpass  Filter  (Figure  5.11) 
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C . Results  of  the  Noise  Analysis 

c2 

The  total  noise  variance  _J1  (equation  (5.7))  was  determined  for 

°0 

each  of  the  filter  realizations  described  in  Section  5.B.  Also,  the 

normalized  power  in  the  passband  of  each  filter  was  computed  by 

a2 

(5.9)  where  ^ anfl  were  chosen  to  be  the  upper  and  lower  0.5  dB 
cuttoff  frequencies.  The  results  are  listed  in  Table  5.9  where  the 
total  normalized  noise  variance  and  the  normalized  passband  noise  power 
are  expressed  in  decibels.  The  entries  for  the  MFB  LP80  filter  are 
shaded  to  indicate  that  this  filter  is  not  physically  realizable  since 
it  contains  delay  free  loops. 

Except  for  the  nonrealizable  case,  the  multiple  feedback  structure 
gave  the  best  noise  performance.  For  the  fifth  order  LP8  transfer 
function,  the  total  noise  variance  of  the  MFB  structure  was  10.4  dB 
better  than  the  cascade  form.  For  the  tenth  order  BP1  and  BP2  filters 


it  was  4.4  and  11.2  dB  better  respectively.  It  was  expected  that  for 
the  tenth  order  bandpass  filter  BP2 , the  MFB  structure  would  show  even 


greater  improvement  over  the  cascade  form  than  for  the  fifth  order  LP8 
filter  since  the  BP2  filter  has  10  poles  grouped  fairly  close  together 
near  z=1  while  the  LP8  filter  has  only  5.  The  reason  that  it  actually 
showed  only  0.8  dB  improvement ( from  10.4  to  11.2  dB)  is  believed  to  be 
due  to  the  fact  that  the  transformation  from  lowpass  to  bandpass 
introduced  four  extra  feedback  paths  in  the  filter  which  added  more 


noise  sources  to  the  structure.  The  fact  that  the  MFB  structure  did  not 
show  much  improvement  over  the  cascade  structure  in  the  BP1  case  was 
expected  since  in  Chapter  4 it  was  seen  that  the  direct  form  second 


order  section  has  the  least  noise  variance  when  its  poles  are  at  +90c 
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Coupled 

Coupled 

Filter 

Cascade 

MFB 

Form  A 

Form  B 

o2/o2 
Lp3  n 0 

26.9 

3^.7 

16.5 

23.5 

19.2 

26.4 

19.3 

26.5 

b 2 

49.1 

y 2 3.9  7 

22.7 

23.0 

LP80  _ ° 

■/ 

°\x 

73.8 

/;44.7 

V/ 

45.4 

45.7 

o2/o2 
n 0 

23.9 

19.5 

27.8 

BP1  „ ^ 

2 .2 
PB/  0 

40.0 

26.3 

34.8 

w 

BP2  _n  ® 

43.0 

31.8 

37.0 

a2  /a2 

PB  0 

50.3 

38.6 

44.6 

Table  5.9. 


Normalized  Noise  Variance 


°n 

3 


and 


Normalized  Passband  Noise  Power 


in  dB. 


Filter 

Cascade 

MFB 

Coupled 
Form  A 

Coupled 
Form  B 

LP8 

5 

6 

5 

7 

LP80 

5 

10 

5 

7 

BP1 

10 

15 

10 

BP2 

10 

15 

10 

Table  5.10.  Number  of  Full  Precision  Multipliers 
Required  For  Each  Filter. 


* 
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It  is  interesting  to  note  however  that  tne  normalized  passband  noise 
power  of  tne  cascaded  second  order  sections  for  BP1  is  13.7  dE  worse 
than  that  of  the  MFB  filter.  This  indicates  that  while  the  total  noise 
variance  of  the  cascade  form  is  not  much  worse  than  that  of  the  MFB  form 
for  b PI,  much  more  of  its  nc ise  power  is  concentrated  in  the  passband . 

The  coupled  forms  also  performed  well  as  expected  from  Chapter  4 
for  the  lowpass  filters.  Coupled  form  A was  7.7  dB  better  tnan  the 
cascade  form  for  LPd  and  26.4  dB  petter  for  the  very  narrowband  LPdO 
filter.  Coupled  form  B was  slightly  worse  than  form  A as  expected 
making  form  A preferable  to  form  B since  it  contains  fewer  multipliers. 
Tne  MFB  filter  was  2.7  dB  better  than  coupled  form  A for  LPd  wnich  is 
reasonable  since  it  is  observed  that  the  network  structure  of  tne  LPd 
coupled  form  filter(Figure  5.3)  is  exactly  that  of  the  MFB  filter(Figure 
5.6)  with  every  other  feedoack  path  eliminated.  The  extra  feedback 
paths  in  the  MFB  filter  allow  more  free  parameters  in  its  design.  The 
performance  of  coupled  form  A for  the  BP2  filter  fell  halfway  between 
that  of  the  MFB  structure  and  that  of  the  cascade  structure  showing  that 
for  poles  near  z=1  , this  form  i3  superior  to  the  cascade  filter.  For 
the  BP1  filter,  coupled  form  A was  3.9  dB  worse  than  the  cascade  form 
which  was  expected  since  the  coupled  form  did  not  perform  as  well  as  the 
direct  form  for  poles  near  the  imaginary  axis  in  Chapter  4.  However  the 
normalized  passband  noise  power  of  the  cascade  filter  was  5.2  dB  worse 
than  that  of  the  coupled  form  filter,  again  indicating  that  tne  cascade 
form  has  much  more  noise  power  concentrated  in  the  passband.  Finally  it 
is  seen  that  the  nonrealizable  LPfiO  MFB  filter  had  about  the  same  noise 
variance  as  the  coupled  form  although  the  result  of  the  LP3  filter 
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design  indicates  that  better  noise  performance  can  be  obtained.  This 
discrepancy  is  thought  to  lie  in  the  fact  that  the  frequency 
transformation  used  to  obtain  the  filter  introduced  four  extra 
multipliers  which  count  as  noise  sources.  This  filter  was  included  for 
comparison  only  and  would  have  to  be  redesigned  to  eliminate  the  delay 
free  loops  for  a more  meaningful  analysis.  Such  a design  could  result 
in  lower  output  noise. 

JL  liardwgrg  Considerations 

Table  5.10  compares  the  number  of  full  precision  fixed  point 
multipliers  required  for  the  filters  of  Table  5.9.  Since  the  cost  of 
these  multipliers  is  a major  portion  of  the  cost  of  a filter  it  is 
desired  to  use  as  few  as  possible  while  at  the  same  time  reducing  the 
noise  power  output  of  a filter.  The  noise  output  of  a filter  decreases 
by  about  6 dB  for  each  bit  added  to  the  register  length  of  a filter[6]. 
Therefore  for  a given  signal-to-noise  ratio,  filter  structures  with 
lower  output  noise  can  be  realized  with  shorter  register  lengths.  From 
Table  5*9,  it  can  be  seen  that  the  lowpass  LP8  MFB  filter  is  better  than 

the  coupled  form  filters  by  le-..,  than  1 bit.  However  the  MFB  filter 

requires  one  more  precision  multiplied 3 multipliers  are  unity)  than 
coupled  form  A for  a fifth  ordr-  filter.  The  choice  between  the  two 
structures  here  is  not  clear  cut  however  both  offer  between  a one  and  a 
two  bit  improvement  over  the  cascade  structure  for  the  LP8  case  and  a 

five  bit  improvement  for  the  LP80  case.  If  the  LP80  MFB  filter  is 

redesigned  to  eliminate  the  extra  multipliers,  it  may  be  possible  to  get 
a six  bit  Improvement  here.  For  the  BP1  filters,  the  MFB  structure 


offers  less  than  a one  bit  improvement  over  the  cascade  form  while 
requiring  five  more  full  precision  multipliers.  In  this  case  it  might 
be  better  to  add  one  bit  to  the  cascade  structure  rather  than  to  use  the 
MFB  structure.  For  the  BP2  filter  the  MFB  structure  is  two  bits  better 
than  the  cascade  form  and  the  decision  is  not  clear.  It  should  be  noted 
that  the  bandpass  MFB  filter  design  method  used  by  Long  would  contain 
only  one  feedback  multiplier  for  each  stage  instead  of  two  as  in  Figure 
3.3  however  it  is  not  known  if  this  method  will  result  in  similar 
performance.  Coupled  form  A is  one  bit  better  than  the  cascade  form  for 
the  BP2  filter  and  might  be  used  since  it  does  not  require  any  more  full 
precision  multipliers  than  the  cascade  form. 

Finally,  it  was  found  that  use  of  the  scaling  technique  of  Section 
2.C  resulted  in  a savings  in  noise  power  of  from  one  to  three  dB  for  the 
filters  constructed.  Thus  unless  noise  reduction  is  a much  greater 
factor  than  cost,  the  technique  would  not  be  very  useful  here. 
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6.  CONCLUSIONS 

The  transformed  multiple  feedback  structures  performed  well  with 
the  proven  LP8  MFB  filter  structure  used  as  a prototype.  The  best 
performance  relative  to  the  cascade  structure  was  achieved  when  the 
poles  were  grouped  near  z=1  as  in  the  LP8,  LP80,  and  BP2  filters.  The 
MFB  structure  was  also  better  than  the  cascade  form  for  the  BP1  filter 
whose  poles  are  near  z=+j  although  the  improvement  was  not  as  great. 
The  problem  with  using  the  frequency  transformations  of  Section  3.C  on 
the  MFB  filter  is  that  the  transformed  filters  are  physically  more 
complex  than  their  prototypes.  Also,  if  the  transformation  attempts  to 
change  the  filter  bandwidth,  feedback  paths  are  introduced  which  make 
the  filter  physically  unrealizable.  Thus  the  design  of  MFB  filters 
cannot  be  accomplished  in  general  by  these  frequency  transformation 
techniques.  The  results  of  Chapter  5 indicate  however  that  MFB  filters 
can  be  designed  to  outperform  corresponding  cascade  filters(direct  form 
2nd  order  sections)  for  any  pole  locations.  More  research  should  be 
concentrated  on  finding  a design  technique  to  minimize  the  noi3e  of 
multiple  feedback  filters. 

The  coupled  form  second  order  section  was  shown  in  Chapter  4 to 
have  less  noise  variance  than  the  direct  form  for  poles  near  z=+l.  In 
Chapter  5,  filters  constructed  with  cascaded  coupled  form  sections  also 
performed  well  for  poles  in  this  region  when  compared  to  fil  i-ers 
constructed  with  cascaded  direct  form  sections.  For  narrowband  lowpass 
filters(LP8  and  LP80)  the  performance  of  the  coupled  form  approached 
that  of  the  MFB  structure.  It  is  interesting  to  note  that  the  modified 


I 

lowpass  MFB  structure  reduces  to  the  cascaded  coupled  form  structure 
when  every  other  feedback  path  is  removed.  The  simplicity  of  coupled 
form  A combined  with  its  noise  properties  make  this  an  attractive  filter 
structure  to  use  when  the  filter  poles  are  near  z=±1. 

For  both  the  MFB  and  coupled  form  structures,  the  technique  of 
selecting  certain  free  parameters  to  be  unity  resulted  in  improved  noise 
performance  due  to  the  elimination  of  noise  sources.  It  was  shown  in 
Chapter  4 that  the  performance  of  a filter  designed  this  way  can 
actually  be  better  than  that  of  a filter  whose  noise  variance  has  been 
minimized  with  respect  to  the  free  parameters.  This  method  of  reducing 
roundoff  noise  is  also  appealing  because  the  number  of  costly 
multipliers  required  can  be  reduced.  Where  free  parameters  are 
available  in  a filter  structure,  the  investigation  of  this  technique  for 
noise  reduction  may  be  profitable. 
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